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1.  Introduction 


Significant  technical  strides  have  been  accomplished  in  designing  cm- scale  micro  air  vehicles 
(MAVs).  One  approach  for  accomplishing  this  mission  is  to  develop  a  biologically  inspired 
flapping  wing  insect  that  can  maneuver  into  confined  areas  and  possess  hovering  capabilities. 
Insect-like  MAVs  have  great  potential  for  limited  duration  electronic  surveillance  and  sensor- 
based  detection  for  both  military  and  civilian  missions  (/).  Some  advantages  of  MAVs  include 
portability,  low  noise,  low  altitude  aerial  surveillance,  low  radar  cross  section,  low  power 
consumption,  and  ease  of  generation  of  lift  or  thrust  at  low  weight.  MAVs  are  implemented  as 
fixed,  rotating,  or  flapping  wing  solutions.  Fixed  wing  MAVs  use  the  wing  span,  the  angle  of 
attack  as  a  function  of  forward  airspeed,  and  the  surface  area  to  gain  lift  and  thrust;  meanwhile, 
flapping  wing  systems  generate  lift  through  the  flapping  wing  beat  frequency.  A  significant 
advantage  of  the  flapping  wing  systems  is  lower  bandwidth  where  the  cyclic  control  inputs  to 
maintain  lift  and  stability  operate  at  lower  frequencies  as  compared  with  a  rotor  system  (2). 
However,  due  to  the  size  limitations  of  MAVs  (typically,  less  than  15  cm  in  dimension  and  less 
than  20  g  in  mass)  it  is  a  challenge  for  the  flapping  wing  vehicle  to  perform  in  adverse  weather 
conditions  (i.e.,  wind  gust  and  precipitation).  One  tactic  is  to  provide  adaptability  or  shape 
control  of  the  wing  structure  that  can  create  effective  aerodynamic  forces  (5). 


Analysis  of  insect  flight  indicates  that  in  addition  to  the  bending  excitation  (flapping), 
simultaneous  excitation  of  the  twisting  degree-of-freedom  (pitching)  is  required  to  manipulate 
the  control  surface  adequately  (2).  Traditionally,  bimorph  piezoelectric  Pb(Zr,Ti)03  (also 
referred  to  as  PZT)  actuators  have  been  used  in  many  applications  to  excite  the  bending  degree- 
of-freedom.  A  common  bimorph  configuration,  as  shown  in  figure  1,  consists  of  two  thin 
ceramic  plates  bonded  together  and  driven  with  opposing  electrical  fields.  One  plate  expands 
while  the  other  contracts,  resulting  in  a  bending  deflection  of  the  bimorph  construction.  As  a 
result,  a  flapping  motion  is  generated  from  operating  at  the  fundamental  resonance  frequency  of 
the  system  (4). 


Figure  1 .  Conventional  bimorph  actuator. 


1 


In  prior  work,  a  flapping  mechanism  was  developed  that  couples  a  3  cm  PZT-5H  bimorph 
actuator  oscillating  at  resonance  frequency  with  flexural  wing  joints  placed  at  the  distal  end  of 
the  bimorph  (see  figure  2).  However,  in  the  current  approach  bend-twist  coupling  arises  from 
anisotropic  material  symmetry.  In  laminated  or  layered  structures,  bend-twist  coupling  is 
governed  by  the  existence  of  at  least  one  anisotropic  layer  not  aligned  with  the  primary  plate 
axes  (5).  By  adding  to  a  bimorph  PZT  actuator  a  layer  of  off-axis  PZT  segments  active,  thereby 
producing  a  layered  structure  to  be  referred  to  as  a  functionally-modified  bimorph,  bend-twist 
coupling  may  be  introduced  to  the  flexural  response  of  the  layered  PZT.  Furthermore,  by 
selectively  charging  off-axis  layers  in  specific  combinations  with  the  bimorph,  the  response  of 
the  functionally-modified  bimorph  may  be  tailored  yielding  a  biaxial  actuator  to  actively  control 
the  flapping  wing  response. 

PZT  Bimorph  Actuator 


Figure  2.  Flapping  wing  device. 

Several  studies  have  considered  the  use  of  piezoelectric  actuators  to  drive  the  bending  degree-of- 
freedom  of  cantilevered  beams  for  applications  ranging  from  control  surfaces  of  aerospace 
structures  (<5)  to  smart  rotorcraft  blades  {7,8).  Piezoelectric  cantilevered  bimorph  actuators  for 
optical  deflection  and  modulation  in  MEMS  applications  have  been  investigated  by  applying  an 
Euler-Bemoulli  beam  analysis  to  clamped  bender  configurations  (9).  Studies  have  also  been 
presented  that  use  energy  methods  {10)  and  elasticity  solutions  (11)  to  develop  quasi-static 
bending  solutions  that  yielded  parametric  relationships  between  tip  deflection,  applied  loads  and 
moments,  actuator  geometrical  dimensions  and  material  properties.  The  present  study  develops 
a  numerical  solution  to  the  three-dimensional  bending  analysis  of  a  cantilevered  functionally- 
modified  bimorph  to  investigate  the  effects  of  angled  segment  geometry  and  orientation  on  the 
bimorph  bending  response. 
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2.  Geometry  and  Nomenclature 


A  conventional  PZT  bimorph  cantilever  configuration  functionalized  by  the  addition  of  angled 
PZT  segments  shown  in  figure  3  was  considered.  The  coordinate  axis  referred  to  as  x-y-z  is 
shown  in  figure  4,  where  the  origin  is  fixed  at  the  clamped  end  of  the  functionalized  bimorph  at 
the  base  of  the  bottom  layer.  For  the  purpose  of  analysis,  the  functionally- modified  bimorph  was 
idealized  as  a  homogeneous  base  with  homogeneous  angle  segments  attached  on  top,  consisting 
of  transversely  isotropic  PZT-5H  material  (Young’s  modulus,  E  =  65  GPa).  The  stainless  steel 
shim  depicted  in  figure  3  was  ignored 


The  distance  from  the  origin  to  the  midplane  of  the  base  along  the  z-axis  was  h.  The  total 
thickness  of  the  base  bimorph  was  2 h,  and  the  thickness  of  the  angled  segments  was  h.  The 
width  b  shown  in  figure  3  relates  to  the  thickness  h  as  follows:  b  =  0. 1  h.  The  angle  of  the  PZT 
segments  was  9  as  shown  in  figure  4.  The  width  of  each  angled  segment  measured  along  the 
x-axis  was  d,  as  shown  in  figure  5.  The  overall  length  Lt  of  the  functionalized  bimorph  included 
five  angled  segments  separated  by  distance  s  as  shown  in  figure  4. 


Figure  4.  Coordinate  system. 
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In  order  to  develop  section  properties  and  geometrical  relationships  useful  for  design  of  PZT 
configurations,  a  representative  element  of  the  functionalized  bimorph  was  defined  to  consist  of 
a  single  angled  segment  affixed  to  a  bimorph  layer  shown  in  figure  5. 


C  O  A 


CO  A  Y 


Figure  5.  Representative  element. 


The  representative  element  length,  L,  shown  in  figure  5,  was  divided  into  three  distinct  regions, 
namely,  sections  1,  2,  and  3.  Typical  cross-sections  for  are  shown  in  figure  6,  where  the  shaded 
regions  represent  the  angled  segment.  Using  the  described  geometry  and  the  points  A,  B,  C,  and 
D  and  points  Q,  R,  S,  and  T  shown  in  figure  5,  functional  relations  for  y-axis  locations  along  the 
edges  of  the  angle  segments  were  derived,  and  are  presented  below.  Section  1  (0  <x<d),  was 
bounded  in  y-z  by  the  plane  containing  the  origin  O  and  the  points  A  and  C  at  x  =  0  and  the  plane 
containing  points  Q  and  T  at  x  =  d.  As  x  increased  from  0  to  d,  the  width  of  the  angled  segment 
increased  from  0  to  t  with  respect  to  the  y-axis.  The  section  1  cross-section  typified  by  the  y-z 
plane  at  A  -  A  "in  figure  6  was  characterized  by  the  angled  segment  where  the  y-location  of  edge 
AQ  was  fixed  at  y  =  b! 2  and  interior  edge  AT  was  a  function  of  x: 


Vat(x) 


b  t 

- X 

2  d  . 


(1) 


Section  1:  Plane  A'-A" 

Section  2:  Plane  B'-fl" 

Section  3:  Plane  C’-C" 

0£x<d 

dixiL-d 

L-d<x<L 

Figure  6.  Three  representative  element  typical  cross-sections. 
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Section  2  (d  < x  <  L  -cl)  was  bounded  in  y-z  by  the  plane  containing  points  Q  and  T  at  x  =  d,  and 
at  x  =  L  -  cl  by  the  plane  containing  points  R  and  S.  As  x  increased  from  d  to  L  -  d,  the  width  of 
the  angled  segment  remained  constant  at  t.  The  section  2  cross-section  typified  by  the  y-z  plane 
B  '-B  "  in  figure  6  was  characterized  by  the  angled  segment  where  the  y-locations  of  both  edges 
QR  and  TS  were  functions  of*: 


yQR(x) 


t—b 
L  —  2d 


,  ,  b 

(pc — d)  +  — 


y-is  A)  = 


b—t 
L  -  2d 


(L  -  d  -  x) 


b 

2 


(2) 

(3) 


Section  3  (L  -  d  <  x<L)  was  bounded  in  y-z  by  the  plane  containing  the  points  R  and  S  at 
x  =  L  -  d,  and  the  plane  containing  points  B  and  D  at  *  =  L.  As  *  increased  from  L  -  cl  to  L.  the 
width  of  the  angled  segment  decreased  from  t  to  0  with  respect  to  the  y-axis.  The  section  3 
cross-section  typified  by  the  y-z  plane  C  -  C"in  figure  6  was  characterized  by  the  angled 
segment  where  the  y-location  of  interior  edge  RD  is  a  function  of  *  and  edge  SD  is  fixed  at 
y  =  -bl  2: 

yRD(x)  =  ^(L-x)-^-.  (4) 

a  2 


The  following  relationships  defined  axial  and  transverse  width  parameters  d  and  t, 


d -L-bt&nO 


t  =  tan 


--G 

v2  j 


d. 


(5) 

(6) 


The  geometry  of  the  representative  element  prescribed  constraints  on  the  angled  segment  width 
parameters  d  and  t — namely,  d  >  0  and  t  <  b.  Applying  the  constraints  on  the  width  to  the 
relationships  of  equations  5  and  6  yielded  the  following  bounds  on  0,  thus  defining  angled 
segment  orientation, 

cot~1(2Z?/  L)  <G  <  tanTl(L/  b)  .  (7) 


3.  Analysis 


The  results  of  traditional  bending  analyses  have  been  used  to  relate  bimorph  loading  and 
geometry.  For  the  functionally-modified  bimorph  considered  in  the  present  analysis,  the  cross- 
section  of  the  cantilevered  bimorph  was  not  symmetric  in  respect  to  the  horizontal  (x-y)  plane. 
Therefore,  a  three-dimensional  bending  analysis  was  applied  in  order  to  relate  the  actuator 
loading  and  geometric  properties  to  tip  deflection  (12).  In  order  to  gain  a  complete 
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understanding  of  the  complex  three-dimensional  analysis  the  piezoelectric  effects  will  be  not  be 
considered  here. 


Beginning  with  the  Euler-Bemoulli  assumptions — (1)  the  cross-section  is  infinitely  rigid  in  its 
own  plane;  (2)  the  cross-section  of  a  beam  remains  plane  after  deformation;  and  (3)  the  cross- 
section  remains  normal  to  the  deformed  axis  of  the  beam — it  was  assumed  that  no  torsional  loads 
are  applied  and  that  the  beam  bends  without  twisting.  Finally  making  the  assumption  that  the 
origin  of  the  coordinate  axis  system  was  selected  to  align  with  the  centroid  of  the  beam  as  shown 
in  figure  4,  the  following  set  of  equations  govern  three-dimensional  bending  (12): 


d2 

Kz 

d2v 

+K 

d2w 

d 

r  "| 

dx2 

dx2 

dx2  _ 

=  Py 

+  — 
dx 

\yapx  -d-A 

(8) 

d2 

d2v 

d2w 

d 

r  "i 

dx2 

Wyz 

dx2 

+  HCyy 

dx2 

=  Pz 

+  — 
dx 

\ ZaPx  ~Py\ 

(9) 

where,  v(x)  and  w  (x)  represented  displacements  of  the  beam  cross-section  in  the  y-  and  z- 
directions,  respectively.  Distributed  loads  along  the  x-axis  in  the  x-,  y-,  and  z-directions  were 
defined  as  px(x),  py(x),  and  pz(x),  respectively.  Coordinates  (ya,  z.a)  defined  the  location  in  the  y-z 
plane  of  the  axial  distributed  load  px(x).  Distributed  moments  along  the  x-axis  about  the  y-  and  z- 
axes  were  defined  as  qy(x)  and  q-(x),  respectively.  The  boundary  conditions  for  the  fixed  end  of 
the  cantilevered  bimorph  were  given  by 

v  =  w  =  0;  dv/dx=dw/dx=(y,  (10) 


At  the  free  end,  the  boundary  conditions  were  as  follows  (with  the  shorthand  notation  in  the 
parentheses)  (12): 


Hi^+ h;z  = a  -  y,P,  (K  =  a  -  yf, ) 

Hr.  + K x?  =  -Qy-  Z«F«  (M,  =  Qy  -  V, ) 


d_ 

dx 

d_ 

dx 


Hiil+Kdl« 


dx" 


yz 


dx 


HtP-X  +  K  d2w 


dx" 


yy 


dx" 


Py-baPx-dz]  (Vy  =  P\  ) 

Pz-\?aPx+Cly]  (Vz=Pz) 


(11) 


where  Py  and  Pz  represented  applied  shear  at  the  tip  in  the  respective  y-  and  z-directions.  Applied 
bending  moments  in  the  y-  and  z-directions  were  represented  by  Qy  and  Qz,  respectively.  In  the 
shorthand  notation  (equation  1 1  in  parentheses),  M  and  V  represented  the  internal  moments  and 
shears,  respectively.  The  centroidal  bending  stiffness  coefficients  Hyy,  Hzz,  and  Hyz  were  defined 
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Hyy=Hcyy  +  Sz2c ;  HZ2  =  HCZ,  +  Sy1c\  Hyz=Hcyz  +  Syczc ;  (12) 

where,  7/vv,  H^,  and  //vz  represent  the  bending  stiffness  coefficients  with  respect  to  an  arbitrary 
coordinate  axis: 

Hyy  =  j  Ez2dA;  Hzz  =  j*  Ey2dA;  =  \  EyzdA ;  (13) 

AAA 


where,  E  was  the  Young’s  modulus  and  y  and  z  refer  to  coordinates  of  an  arbitrary  point  in  the 
y-z  plane.  Note  here  that  the  bending  stiffness  coefficients  are  functions  of  beam  cross-sectional 
area  in  the  y-z  plane.  For  the  functionally-modified  bimorph,  as  the  bimorph  geometry  traverses 
the  ,i'-axis,  the  beam  cross-section  in  y-z  is  dependent  on  the  relationships  for  y  as  a  function  of  x 
developed  in  equations  1-4.  Therefore,  the  coefficients  of  the  ordinary  differential  governing 
equations  of  equations  8  and  9  and  the  boundary  conditions  in  equation  1 1  are  functions  of  x,  and 
are  considered  to  be  variable  coefficients.  Generally,  ordinary  differential  equations  with 
variable  coefficients  must  be  solved  numerically. 

The  sectional  stiffness  coefficients  S,  Sy,  and  S-  are  defined  as  follows  (note  that  the  coefficient 
Sy  is  a  function  of  jc): 

5=  J  EdA;  Sy  =  J  EydA;  Sz  =  J  EzdA ;  (14) 

AAA 


the  centroid  of  a  the  beam  cross-section  (yc,zc)  as  shown  in  figure  4  is  defined  in  terms  of  section 
stiffness  coefficients 


Tc 


(15) 


Again,  note  that  the  centroidal  location  yc  varies  with  x-location  along  the  beam.  In  order  to 
perform  the  three-dimensional  bending  analysis  of  the  functionally-modified  bimorph,  the  beam 
was  considered  in  sections  along  its  length  (x-axis)  by  defining  a  representative  element,  as 
shown  in  figure  5.  The  governing  equations  applied  to  the  sections  of  the  representative  element 
shown  in  figure  5,  take  on  the  following  piece- wise  form: 


0 <x  <d : 


dx 1 


d x2 


c  d\  c  d2wx 

:  3  dx  *  dx  , 

Hc  — l-+Hc 

^  dx1  }y  dr 


=  o 

=  0 


(16) 


d <x<L-d : 


dx‘ 


dx' 


Ht 


H 


d\ 

dx2 

d2v2 

dx2 


+  H 


+  H 


d2w2 

K lb? 

d2w~, 


*  dx 2 


=  0 

=  0 


(17) 
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L-d <x<L : 


dx' 


dx~ 


Hi 


d2v 


dx 


2  +  Hr- 


H 


d2v-. 


yz  ,  2 
dx 


+  H 


d2wi 

dx2 

d2w-. 


yy  ,2 

dx 


=  0 


=  0 


(18) 


where,  the  subscripts  on  v(x)  and  w  (x)  represent  the  sections  of  the  representative  element  as 
depicted  in  figure  5.  The  boundary  conditions  for  a  concentrated  load  P  applied  to  the  free  end 
of  the  cantilevered  bimorph  configuration  assumed,  using  the  shorthand  notation  from  equation 
11 


x  =  0  :  vl  =  Wj  =  0;  dvl  /  dx=  dwl  /dx=0; 
x  =  Lt:  My}  =  0 

=  0  (19) 


Additionally,  the  following  continuity  conditions  were  applied  to  allow  the  analysis  to  span  the 
sections  of  the  piece- wise  representative  element,  where  the  displacements  and  slopes  at  the 
junctions  of  the  segments  were  defined  as  being  equal: 

x  =  d  :  \\=v2;  =  w2-  dvy  /dx=  dv2  / dx,  dy\\  !dx=  dn’2  / dx,  (20) 

x  =  L  —  d:  v2=v3;  w2=w3;  dv2/ dx=  dv3  / dx,  dw2  / dx=  dw3  / dx 

In  order,  to  find  the  deflection  of  the  entire  functionally-modified  bimorph,  continuity  conditions 
were  applied  to  a  series  of  representative  elements  with  spaces  of  traditional  bimorph  in-between 
as  shown  in  figure  4.  The  result  of  this  process  was  a  system  of  38  fourth  order  coupled  variable 
coefficient  ordinary  differential  equations  (ODEs).  The  system  of  ODEs  was  solved  numerically 
using  the  built-in  symbol  NDSolve  in  Mathematica  Version  8  (13)  to  find  the  tip  displacement 
in  terms  of  bimorph  geometrical  parameters.  Results  for  the  deflection  of  the  cantilevered 
functionally-modified  bimorph  in  respect  to  angled  segment  geometry  and  orientation  are 
presented  in  the  next  section. 


4.  Results  and  Discussion 


The  angled  segments  of  the  functionally-modified  bimorph  formed  an  angle  0  with  the  clamped 
edge  of  the  bimorph  as  shown  in  the  inset  in  figure  7.  Equations  5  and  6  provided  the 
relationships  between  the  angle  9  and  geometric  parameters  that  govern  length  L,  angled 
segment  widths  b  and  t.  Figure  7  shows  the  upper  and  lower  bounds  on  9  as  defined  by 
equations  5  and  6,  respectively.  The  points  labeled  in  figure  7  as  I- VII  correspond  to  the  case 
numbers  in  table  1. 


* 

Mathematica  is  a  trademark  of  Wolfram  Research,  Inc. 
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Figure  7.  Upper  and  lower  bounds  on  angled  segment  orientation  angle  6. 

Presented  in  table  1  are  three-dimensional  bending  results  for  the  bimorph  configuration  of 
figure  4,  where  b  was  fixed  at  6  mm  and  exposed  to  an  applied  tip  shear  P  =  1.0  N  in  the 
negative  z  (downward)  direction.  Cases  I  and  II  (labeled  in  figure  6)  shared  the  same  6 
orientation  for  different  ratios  of  length  L  and  width  b.  Because  b  was  fixed,  decreasing  L/b  ratio 
indicated  increasing  L;  therefore,  the  length  of  the  representative  element  in  Case  II  was  greater. 
The  maximum  section  displacement  results  occurred  at  the  x  =  Lj  tip  location.  The  results  for 
Cases  I  and  II  predicted  that  maximum  deflection  increases  as  length  L  increases,  indicating  a 
strong  dependence  on  overall  length  Lj. 

Table  1.  Representative  element  results  for  b  =  6  mm  exposed  to  P  =  1.0  N  concentrated 
load  in  negative  z-direction  (downward). 


Case 

L/b 

0 

d/L 

Vmax(KT6m) 

W max  (  1 0  3  III) 

I 

0.9 

34° 

0.35 

-0.355 

-0.0572 

II 

1.3 

34° 

0.48 

-0.161 

-0.121 

III 

1.3 

45° 

0.23 

-0.134 

-0.148 

IV 

1.3 

51° 

0.050 

-0.997 

-0.134 

V 

1.9 

45° 

0.47 

-0.643 

-0.374 

VI 

1.9 

51° 

0.35 

-2.239 

-0.331 

VII 

1.9 

O 

O 

VO 

0.085 

-2.667 

-0.322 

Comparison  of  Cases  III  and  V  of  table  1,  also  labeled  in  figure  7,  further  evidenced  the 
dependence  of  the  maximum  deflection  value  on  bimorph  length.  Cases  III  and  V  shared  the 
same  value  of  6.  The  length  of  the  Case  V  configuration  was  greater  than  Case  III.  The 
maximum  deflection  increased  as  length  increased.  Also,  considering  Cases  IV  and  VI,  where  6 
values  were  similar  between  the  two  cases.  The  length  of  Case  VI  was  greater,  and  maximum 
deflection  for  Case  VI  was  greater.  The  results  in  table  1  clearly  establish  that  there  is  a  strong 
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dependence  of  maximum  deflection  on  overall  bimorph  length.  Namely,  as  length  increases, 
maximum  deflection  increases. 

Figure  8  presents  the  deflection  curves  for  Cases  II,  III,  and  IV,  where  the  9  values  are  given  by 
34°,  45°,  and  51°  respectively.  Each  case,  as  shown  in  table  1,  had  a  L/b  ratio  of  1.3,  implying  an 
overall  length  Lj  =  28  mm.  Here,  as  9  increased,  cl/L  decreased,  indicating  a  decrease  in  the 
width  of  the  angled  segment.  Intuitively,  the  increase  in  9  would  seem  to  imply  an  increase  in 
the  stiffening  capability  of  the  segment.  However,  as  shown  in  figure  8,  Case  III  at  <9  =  45°, 
cl/L  =  0.23,  displayed  the  greatest  deflection  and  was  therefore  the  least  stiff  configuration,  while 
for  Cases  II  and  IV,  both  configurations  displayed  greater  stiffness,  compared  with  Case  III,  as 
segment  width  increased  and  decreased,  respectively. 


Figure  9  shows  the  relationship  of  deflection  and  angled  segment  orientation  0  for  L/b  =  1.3. 

The  data  showed  that  as  #  increased  between  the  bounds  prescribed  by  the  representative 
element  geometry,  the  deflection  reached  a  maximum  at  around  9=  45°  and  decreased  until  it 
reached  the  upper  bound  on  9.  As  0  increased,  the  data  in  figure  9  show  that  the  width  parameter 
cl/L  indicated  decreasing  width  cl.  The  change  in  the  width  of  the  angled  segment  resulted  in  a 
change  in  the  section  bending  stiffness  coefficients  Hyy,  Hzz,  and  //^defined  in  equation  12.  As  9 
increased  to  9=  45°,  the  section  bending  stiffness  coefficients  decreased  to  a  value  that  nearly 
coincided  with  the  section  properties  for  the  traditional  bimorph  (without  the  angled  segments). 
Therefore,  due  to  the  change  in  angled  segment  width  brought  about  by  the  change  in  angle  9 \ 
the  section  properties  influenced  the  bending  response  so  that  the  resistance  to  bending  in  the 
bimorph  decreased.  Thus,  the  orientation  and  width  of  the  angled  segments  interacted  in  such  a 
way  that  the  bending  capacity  of  the  traditional  bimorph  served  as  a  lower  bound  on  deflection. 
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The  increase  in  deflection  after  6=  45°  may  have  been  due  to  error  typical  to  numerical 
solutions. 


Figure  9.  Effect  of  orientation  angle  6  on  functionally-modified  bimorph  deflection  for  L/b  =  1.3. 


5.  Conclusions 


A  quasi-static  three-dimensional  bending  analysis  of  a  functionally-modified  bimorph 
configuration  has  been  presented.  The  analysis  was  conducted  on  an  idealized  configuration 
where  angled  PZT  segments  were  assumed  to  be  affixed  to  a  base  traditional  bimorph 
configuration.  The  presence  of  the  metal  shim  was  neglected  analytically.  In  general,  bimorph 
deflection  increased  with  increase  in  bimorph  length.  For  fixed  bimorph  length  LT,  as  6 
increases,  deflection  increases,  due  to  changes  in  the  bimorph  sectional  stiffness  properties.  The 
prescribed  geometry  of  the  functionally-modified  bimorph  configuration  provided  constraints  on 
the  orientation  angle  of  the  PZT  segments  affixed  to  the  bimorph  surface.  The  orientation  and 
width  of  the  angled  segments  interacted  in  such  a  way  that  the  bending  properties  of  the 
traditional  bimorph  act  as  a  lower  bound  on  deflection.  Future  work  will  include  numerical 
validation  of  the  three-dimensional  bending  analysis  via  finite-element  analysis  and  experimental 
validation,  with  the  ultimate  goal  of  developing  a  design  tool  for  biaxial  actuators  for  cm- scale 
flapping  wing  MAVs. 
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